Critical Thickness in Dewetting Films 
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We study dewetting of thin polymer films with built-in topographical fluctuations produced by 
rubbing the film surface with a rayon cloth. By varying the density of imposed surface defects, 
we unambiguously distinguish spinodal dewetting, which dominates in liquid films thinner than a 
characteristic thickness = 13.3 nm, from heterogeneous nucleation in the thicker films. Invariance 
of this characteristic thickness upon more than a decade change in the defect density makes kinetic 
effect an unseemly origin. A crossover of the spinodal line provides a consistent picture. This 
interpretation, however, contends the current understanding of molecular interactions in apolar 
liquid films. 
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A physical state is unstable if the curvature of its free 
energy, G with respect to the order parameter, h is nega- 
tive. According to the mean-field theory of Cahn , such 
states may proceed to equilibrium via a spinodal mech- 
anism that results in a spatial modulation of the order 
parameter. It can be shown that the mode with wave- 
length, Xm = (-47r27/d2G'/d/l2)5 

grow the fastest, 

where 7 is the rigidity of the system order parameter 
against deviations from homogeneity. Alternatively, the 
system may evolve to equilibrium by nucleation of hetero- 
geneity at defect sites from within. For a long time, the 
distinction between these two mechanisms in the dewet- 
ting of liquid films from a non- wetting substrate has been 
the subject of many contentious debates [p|-pl]|. 

In supported apolar liquid films, the non-retarded 
Lifsitz-van der Waals interactions render the unit area 
free energy the form —A/127Th'^ where A is the 

Hamaker constant and h - the system order parameter 
- is the liquid film thickness. Therefore, thin liquid films 
wherein A > will rupture spontaneously by spinodal 
growth of surface capillary waves, for which the fastest 
growing mode has the wavelength. Am = (1677^7/^) 2 /i^ 
|0,|3|jlJ] (Here, 7 is the surface tension of the liquid.) 
Previous experiments find that dewetting of liquid films 
proceed according to a generic sequence (^-|^ beginning 
with a thickness modulation, often noted as holes, fol- 
lowed by the coarsening of the holes and formation of a 
polygon network, with the eventual breakup of the poly- 
gon network into droplets due to the Rayleigh instability 
. While the latter processes are governed by hydrody- 
namics of dewetting and well understood p^-p^, opin- 
ions are divided on the mechanism (i.e. spinodal or het- 
erogeneous nucleation) underlying the crucial initial step. 
An /i~2 dependence frequently found in the dewetting 
wavevector [§|-0,|o) has been taken as evidence for spin- 
odal dewetting. However, Jacobs et al. |^ found the holes 
to have no positional correlation, thereby unlikely a re- 
sult of spinodal surface undulations. In addition, the 
dependence apparent in the experimental characteristic 
wavevector has also been argued to resemble an exponen- 
tial variation ^ exp(—h/L) [L = constant), fitting to be 
considered an evidence of heterogeneous nucleation from 
intrinsic defects. Furthermore, it was found in dewet- 
ting experiments of metal thin films [p^ and evaporating 
volatile liquid films that morphologies due to both 
kinds of mechanisms may sometimes coexist. This adds 
further complications to the attempt to make a clean dis- 
tinction between the two destabilizing mechanisms using 
conventional approaches that rely on direct analysis of 
the dewetting patterns. 

In this experiment, we used polystyrene (PS) spin- 
coated on silicon covered with an oxide layer as the 
model system, which is the most studied for this problem. 
We artificially introduce small height fluctuations, 5h 
(< 1 nm) to the thin film samples by mechanically rub- 
bing the sample surface with a rayon cloth before dewet- 



ting initiates. For liquid films that dewet by nucleation, 
the characteristic length of the dewetting pattern depend 
on the density of surface fluctuations introduced. But for 
those fllms that dewet by a spinodal process, the char- 
acteristic length should remain unaffected |17|. There- 



fore, simply by comparing the dewetting morphologies 
with and without the rubbing-induced surface defects, 
one can unambiguously identify which of the two dewet- 
ting mechanisms dominates. Our results show that the 
final pattern of all dewet films is always dictated by just 
one mechanism. For the present system (PS/Si02/Si), it 
is found that spinodal dewetting dominates in the thin- 
ner films {h < 13.3 nm) whereas heterogeneous nucle- 
ation dominates in the thicker ones. Abruptness of the 
transition ai h — 13.3 nm between the two mechanisms 
suggests that it is of a thermodynamic rather than ki- 
netic nature. We attribute the transition to a crossover 
of the spinodal line, which however may present chal- 
lenge to our understanding of molecular interactions for 
the present liquid film system. We will elaborate further 
on this point below. By examining the dewetting mor- 
phologies in different regions with respect to the critical 
thickness, we are able to make direct comparison between 
experiment and theoretical predictions on dewetting 
pattern formation. Excellent agreement between theory 
and experiment is obtained. Our results thus resolve the 
controversy by showing that the two mechanisms are op- 
erative in different thickness regions separated by a tran- 
sition. 

The PS sample was purchased from Scientific Polymer 
Products (Ontario, NY). It has a molecular weight of 
13. 7K Da and polydispersity 1.1. The glass transition 
temperature was measured to be ~ 99 °C by differen- 
tial scanning calorimetry. The substrate overlayer was 
prepared by wet oxidation: 4" diameter Si(lOO) wafers 
cleaned as described previously jl^ were annealed at 
1000 °C in a reaction chamber filled with water vapor 
(produced by a simultaneous injection of 5.6 L/min. H2 
and 4 L/min. O2 into the chamber) for 9 mins. 40 sec. 
The thermal oxide was confirmed to be 106 nm thick 
by spectroscopic ellipsometry with better than 1 % uni- 
formity. The oxidized wafers were then cut into square 
pieces of 1 x 1 cm^, and cleaned |Q before coated with 
the polymer. Both the substrate cleaning and polymer 
coating were carried out in a clean room. Before use, 
the polymer films were annealed at 100 °C under vacuum 
(10~2 torr) for 5 h to remove the residue solvent. No sign 
of dewetting was detectable in the samples after anneal- 
ing. To produce the artificial topographical ffuctuations, 
a piece of rayon cloth loaded with a normal pressure of 
10 g/cm^ was rubbed against the film surface at a con- 
stant speed = 1 cm/s by an electric motor. The density 
of surface defects was varied by changing the number of 
rubs applied to the films. Sample topography was char- 
acterized by a Seiko Instruments (Chiba, Japan) SPA- 
300HV atomic force microscope (AFM) . Fig. |^ shows the 
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Fourier spectra for the surface topography of five freshly 
rubbed PS films that had been rubbed with different 
number of times {N = 3 to 15). These spectra were ob- 
tained by radial averaging the two-dimensional (2D) fast 
Fourier transformation (FFT) of AFM topographical im- 
ages (5x5 /xm^) obtained from the samples. As seen, all 
spectra look alike. Except for a gentle, but obivious rise 
towards q — and a small peak at g = 120 /im~^, the 
spectra are overall speaking rather flat. From the FFT 
images (for example, left inset of Fig. [^), one can see 
that the spectra are dominated by anisotropic topogra- 
phies parallel to the rubbing direction. We found that the 
film roughness could be fitted to (0.26^ -I- 0.16^ A^) 2 (nm) 
(solid line in the right inset of Fig. |l|). All these char- 
acterizations about the topography of the rubbed films 
suggest that addition of height fluctations in the films 
by individual rubs are independent events. It thus fol- 
lows that the density of rubbing-induced surface defects 
should increase linearly with N. 

Dewetting experiments were carried out in air. The 
films were annealed at 180 °C for 10 to 60 mins. We de- 
duce the characteristic wavevector, q* of the final dewet- 
ting patterns from the peak of their radial averaged 
Fourier spectra. When the initial holes do not coalesce to 
form a network of polygons before break up into droplets, 
which has been found in thinner films where h < 13.3 nm, 
(g*)^ will be ^ Nd, the areal density of the final liquid 
droplets For the thicker ones, (g*)^ is Np, the areal 
density of the polygons [|-§. Fig. I depicts q* vs. h for 
samples rubbed by different number of times {N = to 
10) in a 3D plot. As seen, the data are divided into two 
groups at h — 13.3 nm: Below 13.3 nm, q* is indepen- 
dent of N; but above, it increases with increasing N. 
This result unambiguously evidences that samples with 
h < 13.3 nm rupture by a spinodal mechanism whereas 
those with h > 13.3 nm rupture by nucleation of holes. 
We note that q*{h) of the unrubbed films fits quite well 

to - /i(-2.1±0.05) foj. ^ < 3 j^^j^ j^j^j _ ;j(-2.5±0.1) 

for h > 13.3 nm (solid lines in Fig. ||). The fact that 
the latter is not notably different from the h^^ depen- 
dence could have been easily mistaken as the signature 
of spinodal dewetting. We have model-fit the data to 
an exponential function (data not shown) and obtained 
similar good fits as Jacobs et al. found 1^. This find- 
ing confirms the earlier suggestion by these authors that 
the q* ^ h? scaling characteristic is insufficient for the 
distinction between spinodal and nucleation dewetting. 

We now examine evolution of the dewetting morphol- 
ogy of unrubbed PS films in different thickness regimes 
about the threshold at h — 13.3 nm. Fig. I shows AFM 
topographical images obtained from a 6.8 nm thick PS 
film upon quenching at different annealing times from 
3 sec. to 23 min. at 145 °C. As seen, uniformly sized 
holes first appeared after 10 sec, which enlarged and de- 
veloped into a bicontinuous structure after 7 mins. In 
Fig. ^, we contrast evolution of the rupturing morphology 



at 145 °C of PS films near the threshold: one just below 
at h — 12.8 nm and one just above at h = 13.8 nm. For 
the h = 12.8 nm thick PS film (Fig. ^(a)), non-uniformly 
sized holes were first found scattering randomly across 
the sample after 4 mins. They grew in size, and more 
holes emerged with time. Until ^ 28 mins. the holes 
filled up the entire area of the film whereupon the thin 
ribbons connecting the holes broke up into droplets. The 
h = 13.8 nm thick PS film, though differ in thickness from 
the h — 12.8 nm one by only 1 nm, undertook a dras- 
tically different pathway (Fig. |^(b)). As seen, randomly 
distributed holes also appeared in the first few minutes, 
but their size and number never grew large enough to 
cover up the entire film at similar times it had taken the 
h = 12.8 nm thick film. Instead, only bigger holes with 
a much lower areal coverage could be seen growing and 
emerging, which after a much longer time of ~ 100 mins., 
finally touch each other and coalesced into polygons be- 
fore the thin ribbons connecting the polygons broke up 
into droplets. 

With illustrations of Figs. ^ and ^, we may summarize 
on what kinds of conformation can be formed in the early 
stage of spinodal dewetting of supported liquid films. Ba- 
sically, both a bicontinuous structure (Fig. ^ and distri- 
bution of circular holes (Fig. ^(a)) can form in the initial 
stage of a spinodal process. Our results show that the 
bicontinuous structure is formed only in samples in the 
deep spinodal region whereas circular holes are formed 
in samples closer to the spinodal line, consistent with 
the 3D nonlinear calculations of Sharma et al. It is, 
however, noteworthy that the linear theory ^ will have 
predicted only the formation of a bicontinuous structure. 
Our result, therefore, clearly demonstrates inadequacy 
of the linear theory in accounting for the spinodal pro- 
cess. Another point to note is that the circular holes 
formed generally do not have positional correlations (see 
Fig. §(a)) although the charcteristic q* is still the same 
as the one predicted by the linear theory ||l|,|l3|,0 , which 
is again in consistency with Ref . 1^ . 

It remains to understand what causes the transition be- 
tween different dewetting mechanisms near h — 13.3 nm. 
Similar transition has also been observed in previous ex- 
periments . Meredith et al. attributed the transition 
to a crossover between spinodal and nucleation process 
due to kinetic dominance of the latter over the former at 
large h |l^. We note that the critical thickness remains 
unchanged even when the density of rubbing-induced de- 
fects is varied by more than an order of magnitude. Ki- 
netic effects will require that the critical thickness should 
change by more than a factor ~ 105(= 1.58) |19], which 
is obviously not evident in the data of Fig. 2r There- 
fore, it is most probable that the transition represents a 
crossover of the spinodal line. This interpretation, how- 
ever, necessitates that the free energy has an inflection 
point at /i « 13.3 nm. Using the more precise form of the 
free energy proposed in Ref. to take into account that 
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the present system, is a PS/Si02/Si tri-layer rather than 
a simple bilayer one for which the G ^ A/K^ relation is 
derived, one will obtain a prediction for the position of 
the spinodal line to be at /i > (5 = 106 nm, the thickness 
of Si02 in this experiment. At present, the reason for 
the discrepancy between the predicted and the observed 
transition at h = 13.3 nm is not entirely clear. It is 
though apparent that retardation effects may no longer 
be negligible for the Lifsitz-van der Waals interactions 
across the Si02 layer since 5 « 100 nm is quite large in 
this study [Q. However, retardation effects will more 
likely worsen the discrepancy between the predicted and 
the experimental critical thickness. Studies are currently 
being carried out to understand how the Si02 inter-layer 
thickness affects the position of the critical thickness. It 
is hoped that with the insight thus gained about the form 
of van der Waals interactions in a tri-layer system, sta- 
bility of thin liquid films can be premeditated. 
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63, 061603 (2001), it was found that the rubbing- induced 
strain in PS films is fully relaxed within seconds upon 
annealing above the Tg of the polymer. Therefore, any 
perturbation due to the strain on G is relaxed to zero 
well within the spinodal rupture time, thus cannot play 
a major role in the film rupture process. 

[18] X. P. Wang, X. Xiao, O. K. C. Tsui, Macromolecules, 
34, 4180 (2001). 

[19] The growth rate for the size of nucleated holes ~ 7^^/?? 

(where 9 is the equilibrium contact angle between 
the liquid and the substrate surface) is relatively inde- 
pendent of h, but the amplitude growth rate of surface 
undulations, Rm ~ h^^ ||5[ p^ , p^ decreases rapidly as h 
increases. 

FIG. 1. Fourier spectra of surface topography of PS films 
rubbed with difi^erent number of times. Left inset shows FFT 
image obtained from an AFM topographical image of a PS 
film that was rubbed 10 times. Right inset shows the rough- 
ness of rubbed PS films vs. number of rubs. Filled circles are 
data. Solid line is the model fit (see text). 

FIG. 2. 3D plot of the characteristic wavevector, q vs. h 
for PS films rubbed by different number of times from to 
10. Solid lines are fits to a power law. Thin dotted lines are 
contour lines at /i = 7, 13.3 and 30 nm. 
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FIG. 3. AFM topographical images of a 6.8 nm thick PS 
film taken after quenching at different annealing times from 
10 sec. to 23 mins. at 145° C. 

FIG. 4. Comparison between dewetting morphology of (a) 
12.8 nm and (b) 13.8 nm thick PS films at 145° C. All pictures 
are AFM topographical images except the last two in (b), 
which are optical micrographs. 
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